ISSN 1364-0380 (on line) 1465-3060 (printed) 



21 



Qeometry & Topology 
Volume 3 (1999) 21-66 
Published: 22 March 1999 




Classical 6j— symbols and the tetrahedron 



Justin Roberts 



Department of Mathematics and Statistics 
Edinburgh University, EH3 9JZ, Scotland 

Email: justin@maths.ed.ac.uk 



Abstract 

A classical 6j -symbol is a real number which can be associated to a labelling 
of the six edges of a tetrahedron by irreducible representations of SU (2) . This 
abstract association is traditionally used simply to express the symmetry of 
the 6j-symbol, which is a purely algebraic object; however, it has a deeper 
geometric significance. Ponzano and Regge, expanding on work of Wigner, 
gave a striking (but unproved) asymptotic formula relating the value of the 6j- 
symbol, when the dimensions of the representations are large, to the volume 
of an honest Euclidean tetrahedron whose edge lengths are these dimensions. 
The goal of this paper is to prove and explain this formula by using geometric 
quantization. A surprising spin-off is that a generic Euclidean tetrahedron gives 
rise to a family of twelve scissors-congruent but non-congruent tetrahedra. 
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1 Introduction 

A classical 6j'-symbol is a real number which can be associated to a labelling of 
the six edges of a tetrahedron by irreducible representations of SU (2) , in other 
words by natural numbers. Its definition is roughly as follows. 

Let V a (a = 0, 1, 2, ... ) denote the (a + 1) -dimensional irreducible represen- 
tation. The SU{2) -invariant part of the triple tensor product V a ® Vj, <8> V c is 
non-zero if and only if 

a<b + c b<c + a c<a + b a + b + c is even (1) 

in which case we may pick, almost canonically, a basis vector e abc (details are 
given below). 

Suppose we have a tetrahedron, labelled so that the three labels around each 
face satisfy these conditions: we will call this an admissible labelling. Then we 
may associate to each face an e-tensor, and contract these four tensors together 
to obtain a scalar, the 67 -symbol, denoted by a picture or a bracket symbol as 
in figure 1. 




Figure 1: Pictorial representation 

This tetrahedral picture is traditionally used simply to express the symmetry 
of the 67— symbol, which is naturally invariant under the full tetrahedral group 
S4. However, it has a deeper geometric significance. To an admissibly-labelled 
tetrahedron we may associate a metric tetrahedron r whose side lengths are 
the six numbers a,b, . . . , f . Its individual faces may be realised in Euclidean 
2-space, by the admissibility condition (1). As a whole, r is either Euclidean, 
Minkowskian or flat (in other words has either a non-degenerate isometric em- 
bedding in Euclidean or Minkowskian 3 -space, or has an isometric embedding 
in Euclidean 2-space), according to the sign of a certain polynomial in its 
edge-lengths. If r is Euclidean, let a ,0b, ■ ■ ■ ,0f be its corresponding exterior 
dihedral angles and V be its volume. 
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Theorem 1 (Asymptotic formula) Suppose a tetrahedron is admissibly la- 
belled by the numbers a,b,c,d,e,f. Let k be a natural number. As k — ► oo, 
there is an asymptotic formula 

ka kb kc 



kd ke kf 



cos 



|y~](fca + 1)^ + ^ j if t is Euclidean 



37rVfc 3 

exponentially decaying if r is Minkowskian 

(2) 

(where the sum is over the six edges of the tetrahedron). 

A (slightly different) version of this formula was conjectured in 1968 by the 
physicists Ponzano and Regge, building on heuristic work of Wigner; they pro- 
duced much evidence to support it but did not prove it. It is the purpose of this 
paper to prove the above theorem using geometric quantization, and to explain 
the relation between SU(2) representation theory and the geometry of M 3 . 

The formula has a lovely and peculiar consequence in elementary geometry. It 
is well-known that a generic tetrahedron is not congruent (by an orientation- 
preserving isometry of ]R 3 ) to its mirror- image, but is scissors-congruent to it 
(in other words, the two tetrahedra are finitely equidecomposable). Inspired 
by the additional algebraic Regge symmetry of 6j -symbols and the asymptotic 
formula above, one may derive from a generic tetrahedron a family of twelve 
non-congruent but scissors-congruent tetrahedra! 

Section 2 contains the algebraic and section 3 the differential-geometric pre- 
liminaries. Section 4 is a warm-up example, computing asymptotic rotation 
matrix elements for SU (2) representations. It works in the same way as the 
eventual computation (in section 5) for the 67— symbol, but is much simpler and 
displays the method more clearly. Section 6 contains the geometric corollaries 
mentioned above and further notes on the Ponzano-Regge paper. 

Throughout the paper, the symbol "~" denotes an asymptotic formula, whereas 
denotes merely an approximation. 
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time collaborating with John Barrett, trying to find a good method of doing 
the actual calculations required. Neither of us had much success during this 
period, and the details presented here were worked out by me later. (I feel quite 
embarrassed at ending up the sole author in this way.) I am especially grateful 
to John for many lengthy, interesting and helpful discussions on the subject, 
and also to J0rgen Andersen, Johan Dupont, James Flude, Elmer Rees, Mike 
Singer and Vladimir Turaev for other valuable discussions. 
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2 Definition and interpretation of 6j— symbols 
2.1 Combinatorial definition 

The simplest definition is via Penrose's spin network calculus, which is related 
to Kauffman bracket skein theory at A = ±1. The details are in the book 
of Kauffman and Lins [6] . There is a topological invariant ( ) of planar links 
(systems of generically immersed curves) defined by sending a link L to 

/^\ f ^number of loops(L) / ^number of crossings(L) ^g-j 

It extends to an invariant of suitably-labelled trivalent graphs in S 2 , for example 
the Mercedes (tetrahedron) and theta symbols shown in figure 2. To define it, 




Figure 2: Mercedes and theta graphs 

we replace each edge by a number of parallel strands equal to its label, and 
connect them up without crossings at the vertices (this imposes precisely the 
conditions (1) on the the three incident labels). Then we replace this diagram 
by the set of all planar links obtained by inserting a permutation of the strands 
near the middle of each edge. Finally, evaluate each of these using (3), add up 
their contributions, and divide by the number of such diagrams (the product 
of the factorials of the edge-labels). Explicit evaluations of these quantities are 
given in [6]. 

Definition 2 The 6j -symbol shown in figure 1 is defined as the spin-network 
evaluation of the above admissibly-labelled Mercedes symbol, divided by the 
product of the square-roots of the absolute values of the four theta symbols 
associated with its vertices. It is manifestly S4 -invariant. 

Remark It is important to note that the spin-network picture is dual to the 
one drawn in figure 1. There, the trilinear invariant spaces are associated with 
faces of the tetrahedron, whereas in the Mercedes symbol they are associated 
with vertices. 
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Although this definition is the simplest, we will need a more algebraic version 
where the 6j-symbol is exhibited as a hermitian pairing of two vectors. 

2.2 Algebraic definition 

Let V\ be the fundamental 2-dimensional representation of SU(2), which we 
will consider as the space of linear homogeneous polynomials in coordinate 
functions Z and W . Then the other irreducibles, the symmetric powers V a = 
S a V , a = 0, 1, 2, ... , are the spaces of homogeneous polynomials of degree a. 
The dimension of V a is a + 1 ; when a is even, it is an irreducible representation 
of SO{3). 

Making Z , W orthonormal determines an invariant hermitian inner product 
(— , — ) on Vi, and induces inner products on the higher representations V a , 
thought of as subspaces of the tensor powers of V\ . The fundamental represen- 
tation has an invariant skew tensor Z®W ' — W ®Z , which induces quaternionic 
or real structures on the V a , according as a is odd or even. 

The SU{2) -invariant part of the tensor product of two irreducibles V a © Vj, is 
zero unless a = b, when it is one-dimensional. Similarly, the invariant part of 
the triple tensor product V a © VJ, © V c of irreducibles is either empty or one- 
dimensional, according to the famous conditions (1). (The meaning of the parity 
condition is clear from the fact that the centre of SU (2) is the cyclic group 7* 2 . 
The other conditions, often written more compactly as \a — b\ < c < a + b, are 
more surprising. Why the existence of a Euclidean triangle with the prescribed 
sides should have anything to do with this will be explained shortly.) 

We want to pick well-defined basis vectors e aa and e abc for these spaces of 
bilinear and trilinear invariants. Since each such space has a hermitian form 
and a real structure, we could just pick real unit vectors, but this would still 
leave a sign ambiguity. To fix this we may as well just write down the invariants 
concerned. Consider the vectors corresponding to the polynomials 

{Z{w 2 - w x z 2 ) a {ZrW 2 - w 1 z 2 f{z 1 w z - w 1 z 3 ) j (z 2 w 3 - w 2 z 3 y 

on C 2 © C 2 and C 2 © C 2 © C 2 respectively, where i = (b + c - a)/2, j = 
(a + c — b)/2, k = (a + b — c)/2. The required vectors are obtained from these 
by rescaling using positive real numbers, to obtain e aa with norm y/cT+T and 
e abc with norm 1. 

Definition 3 Given six irreducibles V a , . . . ,Vf, one can form e abc © e cde © 
e e f a (g) e f db ( SU pposing these all exist) inside a 12-fold tensor product of irre- 
ducibles. One may always form e aa © e bb © • • • © , and permute the factors 

Geometry & Topology, Volume 3 (1999) 



26 



Justin Roberts 



(without reversing the order of the paired factors) to match. Then the hermi- 
tian pairing of these two vectors (inside the 12-fold tensor product) defines the 
associated 6j -symbol by 



where ^ a is simply the sum of the six labels. One should think of it as 
a function of six natural numbers a,b, . . . ,/, defined whenever the triples 
(a, b, c), (c, d, e), (e, /, a), (/, d, b) satisfy the triangle and parity conditions (1), 
in other words when the associated tetrahedron labelling is admissible. It is 
standard to extend the definition to all such sextuples by setting the 6j-symbol 
to zero elsewhere. 

Lemma 4 These two definitions agree. 

Proof (Sketch) The Mercedes spin network evaluation used in definition 2 
can be reinterpreted as an explicit tensor contraction, using Penrose's diagram- 
matic tensor calculus (see [6]). The invariant (L) of a planar link may be evalu- 
ated by making the link Morse with respect to the vertical axis in M 2 , replacing 
cups and caps with i times the standard skew tensor (Z ®W — W ® Z 6 C <8> C 
and its dual, crossings with the flip tensor, and composing these morphisms to 
obtain a scalar. If we draw the Mercedes graph as in figure 3 and use this recipe 
to compute it, we see that it is given as the composition of a vector in y® 2 ^ a 
(coming from the cups in the lower half of the diagram), a tensor product of 
twelve Young symmetrisers (coming from the flips associated to the crossings 
which are introduced where the labels are) and a vector in the dual of V® 2 ^ ° 
(coming from the caps). This can be interpreted as a bilinear pairing between 
a vector in the tensor product of twelve irreps and one in its dual, if we use 
the symmetrisers to project to these. Reinterpreting it as a hermitian pairing 
and including the normalisation factors gives the purely algebraic formulation 
of definition 3. □ 

Remark This definition does not depend on the choice of hermitian form, 
coordinates, or real structure. It does depend on the sign conventions, but these 
can be seen to be sensible (in that the resulting 6j-symbol is S4 -invariant) 
using the lemma. 

Remark A third way of defining the 6j-symbols is to build a basis for (V a 
Vb <8> V c <8> Vd) su( - 2 ^ out of the trilinear invariants using an isomorphism such as 



(V a ® V b ® V c ® V d ) su ^ * 0(K <g> V b ® V e ) su ^ ® (V e ® V c V d ) su ^ 




e 
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Figure 3: Writing the 6j -symbol as a pairing 

where e runs through all values such that (a, b, e) and (e, c, d) satisfy (1). There 
are three standard ways of doing this, corresponding to the three pairings of 
the four "things" a,b,c,d, and the change-of-basis matrix elements are (after 
mild renormalisation) the 6j -symbols. Using this definition makes the Elliott- 
Biedenharn identity (pentagon identity) for 6j'-symbols very clear, but disguises 
their tetrahedral symmetry; therefore we will not consider this method here. 
See Varshalovich et al [15] for this approach. Their definition coincides with 
the two given here, and with the one in Ponzano and Regge (though in these 
physics-oriented papers, half- integer spins are used). 

2.3 Heuristic interpretation 

The representation theory of SU (2) is well-known to physicists as the theory of 
quantized angular momentum. The fundamental 2-dimensional complex rep- 
resentation V\ can be viewed as the space of states of spin of a spin-^ particle. 
The other irreducibles, the symmetric powers {V a = S a V, a £ N}, are state 
spaces for particles of higher spin; indeed, physicists label them by their asso- 
ciated spins j = \a. Quantum and classical state-spaces are very different: the 
classical state of a spinning particle is described by an angular momentum vec- 
tor in R 3 , whereas in the quantum theory, one should imagine the state vectors 
as wave-functions on IR 3 , whose pointwise norms give probability distributions 
for the value of the angular momentum vector. However, when the spin is very 
large, the quantum and classical behaviour should begin to correspond. For 
example, the wave-functions representing states of a particle with large spin j 
should be concentrated near the sphere of radius j in M 3 . 
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Many representation-theoretic quantities, most obviously square- norms of ma- 
trix elements of representations, can be interpreted as probability amplitudes 
for quantum-mechanical observations. Wigner [17] explained the 6j-symbol as 
follows. Suppose one has a system of four particles with spins ^o, ^b, ^c, ^d and 
total spin 0. Then the square of the 6j -symbol is essentially the probability, 
given that the first two particles have total spin ^e, that the first and third 
combined have total spin . (Compare with the third definition in the remark 
in subsection 2.2.) He reasoned that for large spins, because of the concentra- 
tion of the wave-functions, one can treat this statement as dealing with addition 
of vectors in R 3 . Suppose one has four vectors of lengths ^a, ^b, ^c, ^d which 
form a closed quadrilateral. Then, given that one diagonal is ^e, what is the 
probability that the other is \j1 His analysis yielded the formula: 

(a b cl 2 ^ J_ 
\d e fj 3ttV 

where V is the volume of the Euclidean tetrahedron whose edge-lengths are 
a,b,. . . ,f, supposing it exists. He emphasised that this is a dishonest approxi- 
mation: the 6j-symbols are wildly oscillatory functions of the dimensions, and 
his formula is only a local average over these oscillations, true in the same sense 
that one might write: 

For 9^0, cos 2 (kO) « ^ as k -»■ oo. 

Ponzano and Regge improved his formula to one very similar to (2), deducing 
the oscillating phase term from clever empirical analyses, and verified that as 
an approximation it is extremely accurate, even for small irreps. 

3 Geometric quantization 

3.1 Borel-Weil Bott 

To rigorize Wigner's arguments, we need a concrete geometric realisation of 
the representations V a . This is provided by the Borel-Weil-Bott theorem (see 
for example Segal [1] or [3]): all finite-dimensional irreducible representations 
of semisimple Lie groups are realised as spaces of holomorphic sections of line 
bundles on compact complex manifolds, on which the groups act equivariantly. 
We only need the simplest case of this, namely that the irrep V a of SU (2) is 
the space of holomorphic sections of the ath tensor power of the hyperplanc 
bundle C on the Riemann sphere P 1 . If one thinks of these as functions on the 
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dual tautological bundle, which is really just C 2 blown up at the origin, they 
can be identified with spaces of homogeneous polynomial functions on C 2 (ie 
in two variables) with the obvious SU(2) action (or possibly the dual of the 
obvious one, depending on quite how carefully you considered what "obvious" 
meant!) 

Tensor products of such irreps are naturally spaces of holomorphic sections of 
the external tensor product of these line bundles over a product of Riemann 
spheres, for example 

V a ®V b ®V c = ^(P 1 x P 1 x P 1 , C a MC b M C c ) 

with the diagonal action of SU(2) on spheres and bundles. 

The calculation we are going to perform is a stationary phase integration, for 
which we need local differential-geometric information about these holomorphic 
sections. We will take the primarily symplectic point of view of Guillemin and 
Sternberg [4], as well as using the main theorem of their paper (see below). 
Other relevant references are McDuff and Salamon [9] for general symplectic 
background and for symplectic reduction, Kirillov [7] for a concise explanation 
of geometric quantisation, and Mumford et al [12] for the wider context of 
geometric invariant theory. 

3.2 Kahler geometry 

Suppose M is a compact Kahler manifold of dimension 2n. Thus, it has a 
complex structure J acting on the real tangent spaces T p M , a symplectic form 
co and a Riemannian metric B (we avoid the symbol g, which will denote a 
group element). The latter are J-invariant and compatible with each other 
according to the equation 

B{X,Y) =oj(X, JY), 

this being a positive-definite inner product. The Liouville volume form f2 = 
uj n jn\ equals the Riemannian volume form. The hermitian metric on T p M 
(thought of as a complex space) is 

H(X, Y) = B(X, Y) - ico(X, Y) 

which is linear in the first factor and antilinear in the second (the convention 
used throughout the paper). 
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3.3 Hamiltonian group action 

Let G be a compact group acting symplectically on M. We assume that the 
action is also Hamiltonian (ie that a moment map exists) and that it preserves 
the Kahler structure. This will certainly be true in the examples we will deal 
with. 

We let g be the Lie algebra of G . Then an element £ € Q defines a Hamiltonian 
vector field X^, a Hamiltonian and thus a moment map \x: M — > q* 

according to the conventions 



I will later abuse notation slightly and write n(X^) instead of //(£) when I 
want to emphasise the association of the moment map with a vector field cor- 
responding to an infinitesimal action of G, rather than with an explicit Lie 
algebra element. 

3.4 Equi variant hermit ian holomorphic line bundle 

If the symplectic form u represents an integral cohomology class then there is 
a unique (we will assume M is simply-connected) hermitian holomorphic line 
bundle over M, with metric (—,—), whose associated compatible connection 
has curvature form F = (—2m)uj (so that [u] is its first Chern class). 

We assume that G acts equivariantly on C , preserving its hermitian form. The 
space of holomorphic sections V = H (M,jC) is finite-dimensional and has a 
natural left G-action defined by 



(gs)(p) = g-s(g 1 p). 
V becomes a unitary representation of G when given the inner product 



The round bracket notation will be used to distinguish the global or algebraic 
hermitian forms from the pointwise form on the line bundle C, which will be 
written with angle brackets. 

The infinitesimal action on sections is given by the formula 



dn{C) = ix^ (=w(X(,-)). 




& = -(exp(#)s(exp(-#)p)) = (-Vx, + 2iri/i(0)a 



(4) 



This is the fundamental "quantization formula" of Kostant et al. 
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Remark One has to be very careful with signs here, especially as there is 
such variation of convention in the literature. This formula is minus the Lie 
derivative Cx^ s , because we are interested in the left action of G , and the Lie 
derivative is defined using the contravariant (right) action of G on sections via 
pullback. There is an identical problem if one looks at the derivative of the left 
action of G on vector fields, one has: 

^(exp(#)„y) = -£ X( Y = -[Xt,Y] 

provided one uses the standard conventions on Lie derivative and bracket: 

[X, Y] = C X Y, [X, Y]f = X(Yf) - Y(Xf) 

For further comments on sign conventions see McDuff— Salamon [9], remark 3.3, 
though note that we do not here adopt their different Lie bracket convention. 
Anyway, it is a good exercise to check that the formula really does define a Lie 
algebra homomorphism: for this one also needs the standard conventions on 
curvature: 

F(X,Y) = [Vx,V Y ]-V [xx] 

and on Poisson bracket: 

{f,g} = -u(X f ,X g ) 

3.5 Complexification 

In the Borel-Weil-Bott setup, one actually has a complex group G c acting on 
C and M. (Of course it does not preserve the hermitian structure on £, but 
its maximal compact part G does.) The action of g c on sections of C is given 
by 

(ir))s = i(rjs) = (-iV ' x n - 2iT[i(ri))s = {-Vjx n - 2iT{i(ri))s 

the last identity coming because s is a holomorphic section, so is covariantly 
constant in the antiholomorphic directions in TM C : 

^x+ijxs = 

The point of this identity is that it gives us information about the derivatives 
of an invariant section in directions orthogonal to the slice ^i _1 (0). 
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3.6 Example 

The (k + 1) -dimensional irrep of SU(2) is obtained by quantising S 2 with a 
round metric and with an equivariant hermitian line bundle £® k of curvature 
kco, where to is the standard form with area 1. We will always view S 2 as being 
the unit sphere in R 3 . In cylindrical coordinates, its unit area form is then 
("Archimedes' theorem") 

uj = -—dO A dz. 

Let £ be an element in the Lie algebra of the circle so that = 1. The 
moment map for the 1-periodic rotation about the z-axis, generated by the 
vector field = 2ird/d6 is fi: S 2 — > R given by fi = \kz. (Here we identify 
the dual of the Lie algebra with R by letting £ be a unit basis vector.) Thus 
the image of the moment map is the interval [— |, |], and in accordance with 
the Duistermaat-Heckman theorem, the length of the interval equals the area 
of the sphere. 

3.7 Kahler quotients 

Given as above a Kahler manifold M and the Kahler, Hamiltonian action of 
a compact group G, we may form the Kahler quotient Mq, which is just an 
enhanced version of the symplectic (Marsden-Weinstein) reduction. 

Let Mo denote the slice /U _1 (0). The G-equivariance of the moment map 
ri: M — > q* (coadjoint action on the right) means that Mq is G-invariant, and 
consequently that lo(X^,Y) = drt(£)(Y) = for any £ G and Y 6 T p Mq. Let 
us suppose that G acts freely on Mo , since this will be enough for our purposes. 
We will use the symbol Qp = {X^(p) : £ G g} to denote the tangent space to the 
G-orbit at p and similarly, igp will denote {JX^(p) : £ G q} . In fact Qp is the 
symplectic complement to T p Mo at p, and iQp is the (Riemannian) orthogonal 
complement to T p M , because B(Y, JX^) = -uj(Y,X^) = for any Y G T p M , 
and the dimensions add up. 

As a manifold, Mq is just just the honest quotient Mq/G. It inherits an in- 
duced symplectic form <jJq whose pullback to Mo is the restriction of that of 
M . Its tangent space T^Mq at a point [p] (the orbit of a point p G M) may be 
identified with its natural horizontal lift, namely the orthogonal complement of 
QP ^ TpMo at any lift p of [p] . This space may also be described as the orthog- 
onal complement of g c p C T p M . As this subspace is complex, T\^Mq inherits 
both a Riemannian metric and a complex structure by restriction. Hitchin 
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proves in [5] that these induced structures make Mq into a Kahler manifold. 
Starting from C over M , we can also construct a hermitian holomorphic line 
bundle Cq over Mq with curvature —2iriijJG ( m particular, the induced sym- 
plectic form is integral), as in [4]. The bundle and connection are such that 
their pullback to Mq agrees with the restriction of C . 

3.8 Reduction commutes with quantization 

Let Q(M) denote the quantization H°(M,C) associated to a Kahler manifold 
with equivariant hermitian line bundle C (which is suppressed in the notation). 
It is a representation of G, so we can consider the space of invariants Q(M) G C 
Q(M). (Whether we use G or G c here is of course irrelevant.) 

The main theorem in [4] is that there is an isomorphism Q(M) G = Q{M G ). 
There is obviously a restriction map from invariant sections over M to sections 
over Mq, more or less by definition of Mq and Cg, so the task is to show 
injectivity and surjectivity. 

A vital ingredient in their proof of surjectivity is fact that norms of invariant 
sections achieve their maxima (in fact decay exponentially outside of) the slice 
Mq. We will rely on this fact too. Suppose s is a holomorphic G-invariant 
section over M, and consider the real function ||s|| 2 on M. It is certainly 
G-invariant, but not G c -invariant. Following [4] we compute the derivative 

(JA^IHI^^vr/i^lHI 2 

by using the quantization formula (4) and the compatibility of hermitian metric 
and connection 

X\\4 2 = (Vxs,s) + (s,V x s). 
Therefore, if j(t) is the flowline starting at p G Mq and generated by JX V , 




s Rm = - 47r M»7)ll«ll?(t) 



and combining with 

^M^) 7 (t) =B(X v ,X v )>0 
we see that indeed the function ||s|| 2 ( t ) has a single maximum at t = 0, ie on 

Mq. 
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3.9 Refinement of the Guillemin— Sternberg theorem 

We need an addition to the "reduction commutes with quantization" theorem. 
Any space Q(M) is in a natural way a Hilbert space, with inner product defined 



where {—,—) is its line bundle's hermitian form. One might imagine that the 
restriction isomorphism Q(M) G = Q(Mq) is an isometry, but in fact it is not. 
However, asymptotically it becomes an isometry if one redefines the measure 
on Mq , as will be shown below. 

First let us refine the observations about maxima of pointwise norms of sections 
given above. We can repeat the argument using the pointwise modulus of 
{si,s 2 ), and establish that its maxima too are on Mq. Also, we can compute 
the second derivatives of the function (si,s 2 ) m the JX V directions, which span 
the orthogonal complement of TMq. Redoing the above calcuation yields, at 
P e M , 

(JX^.JX T] (s 1 ,s 2 )) p = -AirJX^{n{ri){si,S2)) p = -AirB p {X r) ,X^){si,s 2 ) p 

because the moment map is zero at p. Parametrising a regular neighbourhood 
of Mo as M x exp(iU) for U some small disc about the origin in 0, we see that 
to second order the function satisfies 



To understand the asymptotics, we must first understand what is varying! Let 
A; be a natural number. Then one can consider M with the new symplectic 
form kco; its Liouville form scales by k n (recall dimM = 2n), the moment 
map for its G-action scales by k, its Riemannian metric scales by k, and there 
is a new equivariant hermitian holomorphic bundle C® k over M whose Chern 
form is ku> . If s is a G-invariant section of £ then one can consider s k = s® k , 
which is an invariant section of C® k . We will always write the k explicitly 
to indicate the scaling of forms, so that B, uj, £1, X^ and so on retain their 
original definitions. The new pointwise hermitian form satisfies 



Thus, in view of (5), as k — > 00 the invariant section has pointwise norm 
concentrating more and more (like a Gaussian bump function) on the slice Mo . 
This localisation principle rigorises Wigner's ideas and forms the basis for the 
proof of the Ponzano-Regge formula. 



by 




(si, S2>(p,o ~ (si, s 2 ) p e- 2 ^ x ^ for £ G U. 



(5) 



(4,4) = ( Sl ,s 2 ) k . 
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Theorem 5 Let s±,§2 be G -invariant sections of C over M, and s\,S2 the 
induced sections of £q over Mq. Let a: Mq — > R be the function which 
assigns to a point \p] the Riemannian volume of the G -orbit in M represented 
by \p] , and let d = dim(G) . Then, as k — > oo, there is an asymptotic formula 

(~sl4) = [ <sJ,sS>(* B n) ~ f|V /2 / 

Jm \ z / Jm g 

Proof From basic geometric invariant theory [4] we know Mq/G = M ss /G c , 
where M ss = G c Mq is the set of semistable points, an open dense subset of 
M . We view the left hand integral as an integral over M ss and then integrate 
over the fibres of ir: M ss — ► Mq, which are G c -orbits. It helps to think of 
as the Riemannian (coming from B) rather than the symplectic volume form. 

Suppose \p] G Mq and p G Mq. The fibre can be parametrised via the map 
G x iq — > 7T _1 ([p]) given by (g,i£,) i— > exp(in)gp. This is a diffeomorphism 
because of the Cartan decomposition of G c . (Recall we are assuming that G 
acts freely on Mq.) 

Let if) be the pullback to G x ig of the function log(si, S2) , so that log^, S2) 
pulls back to kip. It is invariant in the G directions but has Hessian form in 
the ig directions (at 0) given by 

(i£.ir].kip)o = -ATrkB p (X ( ,X v ). 

The pullback Riemannian metric on ig is given by (3(i£,ir)) = B(X^, X v ) p . 
Consequently the integral of oik^ over ig is asymptotically given by 

<5i, S 2 ) k p [ s e-W&MdvoLp = (h, s 2 ) k p (^) d/2 • 

(The symbol volp denotes the measure induced by the same metric (3 as appears 
in the integrand; changing coordinates to an orthonormal basis, one obtains the 
a standard Gaussian integral, independent of (3.) 

Finally we integrate over the G orbit, picking up the factor vol(Gp) = o~(\p\). 
Substituting into the original left-hand side and separating the powers of k in 
the correct way finishes the proof. □ 

3.10 Example 

Let us return to example 3.6 and check these formulae. If we have area form 
2kuj then sections of £® 2k are identified as homogeneous polynomials in Z, W 
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of degree 2k, and under the circle action there is an invariant section which one 
can write as 

s k : (Z,W) ' t (Z k W k ) G C. 

In the complement of infinity, trivialise C® 2k using the nowhere-vanishing holo- 
morphic section b 2k corresponding to the homogeneous polynomial W 2k . Let 
C = Z/W be the coordinate on this chart. The pointwise norm of b 2k is 
(1 + |C| 2 ) _fc > because a unit element of the tautological bundle above ( is 
(C, l)/y / (l + |C| 2 ), which gets sent to (1 + |C| 2 )~ fc by the section b 2k .) Thus the 
pointwise norm of s k is |C| fc /(l + |C| 2 ) fc - Under stereographic projection from 
the unit sphere in R 3 

= x + iy 
^ ~ 1-z 

we get 1/(1 + |C| 2 ) = (1 - z)/2 and |C| 2 /(1 + |C| 2 ) = (1 + z)/2. Thus the 
pointwise norm-squared of s k is ((1 — z 2 )/4) k . Its global norm-square is 

|, A r - / \^^-\ k 2k—d9f\dz 



4 / 47T 



1 'l-z 2 \ k l 



- :>!■' I I — — ! -f/.: 



4/2 



= 2JfeB(ife + l,ife + l) 

by definition of the beta function B . Evaluating the beta function in terms of 
factorials gives 

Compare this with the computation of the norm on the reduced space, which 
is a single point: one finds the norm-square of s k at this point to be simply its 
value on the equator z = of the sphere, namely A~ k , so the ratio of the two is 
therefore \firk. If one applies theorem 5 with VL = 2uj then one gets the same 
ratio: the scaling factor is -\fkj2 and the length of the equator (the factor a 
in the formula) is 2tt^2/Att, because the usual spherical area form has been 
divided by Att and multiplied by 2 . 



3.11 Orbit volumes 



It is convenient here to note the formula for the volume of the G-orbit with 
respect to some basis (we will eventually have to roll up our sleeves and perform 
explicit calculations) . 
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Lemma 6 If is a basis for and p is an invariant metric on q, then the 
volume of the G -orbit at p G M is 

Proof Pulling back the metric via the diffeomorphism G — ► Gp, g ^ gp gives 
a metric on G of the form 

/3{g*£,g*v) = Bg P {x^x v ) = b p (x^,x v ) 

whose associated volume form differs from that of p by the given factor. □ 

It is also worth giving a formula for the orbit volume when G does not act freely 
on a space. Suppose as in the previous lemma that p is an invariant metric, 
that the stabiliser at p is T , whose Lie algebra is t C g . 

Lemma 7 If {&} is set of vectors in q which, when projected onto into the 
orthogonal complement t 1 C j, forms a basis for that space, 

Proof Repeat the earlier proof with G/T mapping diffeomorphically to the 
orbit, and using the basis {^{\ for the tangent space of G/T . This gives a 
formula like the above except with voIb p {X^ } on top. However, since the 

vectors £j — £j are in t, they map to zero tangent vectors at p, and we can 
simply remove the hats. □ 



3.12 Stationary phase formulae 

The standard stationary phase formula is as follows: on a manifold M 2n with 
volume form 0, , for a smooth real function / with isolated critical points {p} , 
one has 



2vr\ n ^ e »fc/(p). e f sgn(Hess p (/)) 

V / Hess p (/) 



In our computation of we will actually have a complex function tp, so it is 
probably easier to rewrite/generalise the above formula to such a situation as 

(6) 



Jm \ k J p V- Hess pW 
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where the Hessian is now a complex number, and by the square root we mean 
the principal branch (the Hessian must not be real and positive). Of course, 
with tp = if , / real, this reduces to the previous version. 



4 Warm-up example 

In order to demonstrate more clearly the main points of the calculation to come 
in section 5, we will first work out a simpler case. 

Let V2k be an irreducible representation of 50(3), and let 5* be the circle 
subgroup fixing the z-axis in R 3 . With respect to this subgroup, splits 
into one-dimensional weight spaces indexed by even weights from —2k to 2k. 
We may pick unit basis vectors inside these, uniquely up to a sign (by using 
the real structure of the representation). If g G 50(3) is a rotation, we may 
compute the matrix elements of g with respect to such a basis by using the 
hermitian pairing. Most of these depend on the choices of sign, but the diagonal 
elements, those of the form (v,gv), are independent. Below we will compute 
an asymptotic formula for such a matrix element, when v is the zero- weight 
(5 1 -invariant) vector. 

There are in fact explicit formulae for matrix elements given using Jacobi and 
Legendre polynomials, which are well-known in quantum mechanics. One can 
prove the theorem from these more easily, see for example Vilenkin and Klimyk 
[16]. Also, it could be computed explicitly from the example 3.10. But this 
method demonstrates how to do the calculation without such explicit knowledge 
of the sections. 

If v is a weight vector for Sj. then gv is a weight vector for gS],g~ l = Sg Z , 
the subgroup fixing the rotated axis " gz" . The elements (v,gw) can also be 
thought of as elements of the matrix expressing one weight basis in terms of 
the other. 

Theorem 8 As k — > oo there is an asymptotic formula 

<^'>~>/S~{<* +1 >i + i} 

(k) 

where Vq is a unit zero-weight vector in V 2 fc and (5 > the angle through 
which g rotates the z-axis. 
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Proof As in example 3.6, is the space of holomorphic sections of £® 2k 
over S 2 , with SO(3) acting in the obvious way, symplectic form 2kuj, and 
the subgroup S\ acting with moment map kz. The zero-slice is a circle, and 
the reduced space a single point. Therefore there is a one-dimensional space of 
invariant sections, and such a section will have maximal (and constant) modulus 
on the slice z = 0. 

We choose a section s of £® 2 which is S** -invariant and has peak modulus 
1 where z = 0. (The phase doesn't matter, as noted above.) Then s® h is a 

section of £® 2k , invariant under S^, and also with peak modulus 1 at z = 0. 

(k) 

It does not quite represent a choice of Vq because its global norm is not 1 (we 
have fixed it locally instead). 

We must compute an asymptotic expression for 



The denominator is asymptotically \/nk, as we checked in example 3.10. So 
the main work is computing the integral 



where ip = log(s,gs). 

Now gs is an invariant section for Sg Z , whose moment map is simply the "gz 
coordinate", and whose zero-slice u gz = 0" meets z = in two antipodal 
points N and S. (They are both on the axis of g, and N is the one about 
which g is anticlockwise rotation.) Therefore, outside a neighbourhood of these 
two points, the modulus of the integrand is exponentially decaying, and the 
asymptotic contribution to the integral is just from N and S. In fact these two 
points will also turn out to be the critical points of ij), and we will evaluate the 
integral using the standard stationary phase procedure. 

Let us denote by /x, v the moment maps for S\ and Sg Z acting on the sphere 
with symplectic form 2uj . Let X, Y be the generating vector fields correspond- 
ing to these actions, as shown in figure 4. 

The first derivatives of ip can be calculated as follows: start by computing 

X(s,gs) = (V x s,gs) + (s,V x gs). 

The first term can be simplified to 2mfi(s, gs) via the quantization formula (4), 
because s is invariant under the group corresponding to X . The second term 
is quite so easy, but becomes simpler if we write 





X = pY + qJY 
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Y 



X 

4 



JY 



gz = 



z = 



Figure 4: Localisation regions and generating vectors 



for some scalar functions p,q (generically Y,JY span the tangent space), and 
then expand again: 



At either intersection point, the moment maps are both zero, so the whole thing 
vanishes. Similarly, the substitution 



Applying X and Y again to these formulae gives us the second derivatives at 
the critical points. The fact that ultimately we evaluate where the moment 
maps are zero shows we need only worry about the terms arising from the 
Leibniz rule in which the vector field differentiates them. These are evaluated 
using the following evaluations at N : 

Xfi = Xv = 2lj(Y,X) Yfi = 2u(X,Y) Yv = 0. 

In addition, at N we have p = p' = cos/3 and q = q' = — sin/3, by inspection. 
Thus, with respect to the basis {X, Y} , we have the matrix of second derivatives 



To obtain the Hessian, we have to divide the determinant of this matrix by 
(2uj(X, Y)) 2 , to account for the basis {X, Y} not being unimodular with respect 
to the volume form 2lo. Therefore the Hessian at ./V is 



Xijj = 2iti\i — 2'nipv — 2'nqv. 



Y = p'X — q'jX 



gives 



Y ip = 2irip' ii + 2irq' fi — 2mu. 




(2Tri) 2 e i/3 2i sin /3 = -8tt 2 i sin (3e^ . 
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An easy check shows that the complex conjugate occurs at S . 

The 0-order contributions e k ^( N ^ , e k ^^ must be computed. Both are of unit 
modulus, because of our convention on s. Their phases are easy to compute 
using the unitary equivariance of the bundle. Let h be the rotation through tt 
lying in S\ . It exchanges N and S , and preserves the section s . Thus 

eV'(S) = (s,gs)(S) = (h.h- 1 s(hN), gh.h- 1 s(hN)) = (hs{N), ghs{N)) 
= (s(N),h- 1 gh S (N)). 

Now h~ l gh is clockwise rotation through (3 at N, and therefore acts on the 
fibre of £® 2 , which is the tangent space at N , via e~ %13 . Hence 

e i>{s) = e i(3 and similarly e m) = e -v m 



Finally we can put this all together: 



k 



2vr 



-ikfi 



Jk/3 



+ 



k { 87r 2 i sin /3e*^ y/SnHsmpe-^ 
Dividing by the normalising (s,s) = \fnk gives 



, 2 - cos \(2k + l)- + . 
irk sin (3 \ v J 2 4 ' 



Remark The answer consists of a modulus, coming from the modulus of the 
Hessian's determinant and the normalisation factors for the original sections, 
and a phase, coming partially from the phase of the Hessian and partially 
from the global phase shift (0-order terms). The 7r/4 is a standard stationary 
phase term. It will be possible to identify the terms in the 6j-symbol formula 
similarly. 



Remark Other weight vectors behave similarly. If one repeats example 3.10, 
one finds beta integrals with more '(1 — z)'s than '(1 + z)'s or vice versa, and 
the peak of the integrand also shifts to the correct position: a circle of constant 
latitude whose constant z coordinate equals the weight divided by k. Similar 
formulae for general rotation matrix elements may be obtained. 
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5 Asymptotics of 6j— symbols 
5.1 Geometry of the sphere 

The 6j -symbol arises by pairing two SU(2) -invariant vectors in a 12-fold ten- 
sor product of irreducibles, according to definition 3. Let us first consider the 
geometry of single irreducibles. 

Identify the Lie algebra of SO (3) with R 3 by using the standard vector prod- 
uct structure " x " , so that the standard basis vectors ei generate infinitesimal 
rotations v — > x v in space. Define an invariant metric p (the usual scalar 
product ".") by making them orthonormal. Using this scalar product we iden- 
tify M 3 with the coadjoint space also. This metric gives the flows generated by 
the 'ej 's period 2ir, and so gives the circle T in SO(3) generated by any of them 
length 2-7T. We can think of SO (3) as half of a 3-sphere whose great circle has 
length 4ir; such a sphere has radius r = 2, and therefore volume 2n 2 r 3 = 16ir 2 . 
Therefore £0(3) has volume 8tt 2 , and the quotient sphere S0(3)/T (with its 
induced metric) has volume 4-7T . 

The irrep V a is the space of holomorphic sections of the bundle C a on S 2 . To 
write down explicit formulae for the various structures on the sphere it is better 
to view it as the sphere S 2 of radius a in ]R 3 , instead of as the unit sphere. 
The symplectic form with area a is given by 

U X (V,W) = -^^[x.V.w] 

where x is a vector on the sphere, v, w are tangent vectors at x (orthogonal to 
it as vectors in IR 3 ) and the square brackets denote the triple product 

[x.v.w] = x.(v x w). 

The complex structure at x is the standard rotation 

J x (v) = —X X V. 
a 

The final piece of Kahler structure, the Riemannian metric, is then 
B x (v, w) = lo x (v, J x w) = -£^3 (x x v).(x x w) = 4~3'( x2 ( t; - ti; ) ~ (x.v)(x.w)) 

1 

= - — v.w. 

47TO 

This formula agrees with the fact that with the natural induced metric given 
by B x (v,w) = v.w, the area of S 2 is 4ira 2 . 
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The group G = SO (3) acts on the sphere, preserving its Kahler structure. This 
is a Hamiltonian action with moment map being simply inclusion fi a : S 2 ^ . 
If a is even then G also acts on the corresponding hermitian line bundle, but 
if a is odd one gets an action of the double cover SU(2). (Since we need to 
work primarily with the geometric action, it is SO (3) that wins the coveted 
title "G"!) 

Products of Kahler manifolds have sum-of-pullback symplectic forms, and direct 
sums of complex structures. Their Liouville volume forms are wedge-products 
of pullbacks (there are no normalising factorials; this is one reason for the n! in 
the definition of the Liouville form). The moment map for the diagonal action 
of G on Sq x S 2 x Si is therefore [i{x\, X2, X3) = x\ + xi + X3 G M 3 . We know 
from the discussion earlier that the pointwise norm of an invariant section over 
this space will attain its maximum on the set /i = 0, which is in this case the 
"locus of triangles" {xi + X2 + X3 = 0} . G acts freely and transitively on this 
space, except in the exceptional cases when one of a, b, c is the sum of the other 
two. We can safely ignore this case, because it corresponds to a flat tetrahedron 
(about which the main theorem says nothing). 

We need to pick a section s abc corresponding to the invariant vector e abc , whose 
normalisation was explained in section 2. There is a one-dimensional vector 
space of SU(2) -invariant sections of £ m El £® b M C® c over S 2 a x S 2 x S 2 C . 
First we define s abc uniquely up to phase by setting its peak modulus to be 
1. This section does not represent e abc exactly, because we fixed the norm 
locally, instead of globally. However, it is more convenient for calculations, and 
we can renormalise afterwards. Similarly, let s aa be a section over S 2 x S 2 
with peak modulus 1. It will be concentrated near the zeroes of the moment 
map /J:(xi,X2) = xi + X2, namely the anti-diagonal. Note that SO(3) acts 
transitively on the antidiagonal with circle stabilisers everywhere. To represent 
e aa , this section would have to be renormalised so that its global section norm 
was \/a + 1 . 

Fixing the phase of each section is slightly more subtle. The phase of the 'e's 
was fixed by using the spin-network normalisation. The analogue for sections is 
just the same, viewing the Riemann sphere as P 1 with homogeneous coordinates 
Z and W , and defining the sections just as before. The choice will not actually 
matter until subsection 5.9. 

5.2 A 24— dimensional manifold 

Fix 6 natural numbers a, b, ... / satisfying the appropriate admissibility condi- 
tions for existence of the 6j-symbol. 
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We will work on the following Kahler manifold M: 

Ad = S a x *S*^ x S c x S c x *S*^ x *S* e x S e x S j- x S a x S j- x x *S*^ 

This is taken to lie inside (M 3 ) 12 , and a point in it will be written as a vector 
(21,22, • • • ,212). 

There are three useful actions on M. First there is the diagonal action of G 
on all 12 spheres. It has moment map <j>: M — > R 3 given by 

12 

0(xi,x 2 , • • • ,212) = Xl 3 ^- 
1 

Algebraically, this generates the diagonal action of G on the corresponding 
tensor product of 12 irreducible representations 

V a ®V b ®V c ®V c ® ■ ■ ■ ®V b . 

Secondly, one has an action oiG 4 = GxGxGxG, the first copy acting 
diagonally on the first three spheres, the next on the next three, and so on. 
The moment map for this action is //: M — > (M 3 ) 4 , 

(xi,x 2 , ... , X12) !->■ (x 1 + x 2 + 23, 24 + 25 + x 6 ,x 7 + 2 8 + x 9 ,x w + £11 + X12). 

The section = s abc M s cde M s e f a Kl s^ db is a well-defined (given earlier con- 
ventions on s abc ) invariant section for this action with peak modulus 1. 

Thirdly, we let G 6 act on M , each copy acting diagonally on a pair of spheres 
of the same radius. The moment map (which shows precisely how this works) 
is v: M — > (M 3 ) 6 , 

V>(X 1 ,X 2 , ■■■ , 212) = (Xl + 2 9 , 2 2 + Xl2,X 3 + 2 4 , 2 5 + 2n, 2 6 + X 7 , X 8 + X W ) . 

Then s v = s aa Ms bb M- ■ -Ms^ (after a suitable permutation of its tensor factors, 
so that s aa lives over the first and ninth spheres, for example) is an invariant 
section for this action with peak modulus 1. 

5.3 Proof of theorem 1 

Recall from 3 the definition of the 6j -symbol as a hermitian pairing of two 
vectors, and their normalisations. The corresponding geometric formula, in 
terms of the sections S M , s v just defined is 
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where the product on the right denotes simply (ka+l)(kb+l) ■ ■ ■ (kf+l) .) Our 
convention on peak modulus 1 and phase of the sections mean that (s abc )® k = 
g ka,kb,kc an( j s i m ilarly ^ g aa^k _ g ka,ka ^ so this formula is just the pairing of 
tensor products of these sections, with corrections for the global norms of 
and s u as explained at the end of subsection 5.1. 

To extract the asymptotic formula for the 6j -symbol, we therefore need asymp- 
totic formulae as k — > oo for the three integrals: 

i = (s*,sj) = / <s*,s*>fc 12 n 

JM 
J M 

h = («)= / (s k v ,s$)k 12 n 

JM 

(Note the explicit inclusion of all factors of k; everything else is unsealed.) 
We can in fact immediately write down asymptotic formulae for the correction 
integrals I M , I v using theorem 5, because the reductions M G i , M G <> are both 
single point spaces, and the sections have modulus 1 over these points. 

h - (^vol^-^O)) (8) 

h - (^vol^-^O)) (9) 

(In the second case one must actually reconsider the proof of the theorem, 
because G 6 does not act freely on the set z/ -1 (0), but there is no problem.) 

The remaining integral / is evaluated by reduction to an integral over Mq 
followed by the method of stationary phase, which fills the rest of this section. 

5.4 Localisation of the integral / 

As k — ► oo, the integrand decays exponentially outside the region where both 
moment maps fj,, v are zero, because it is dominated by the pointwise norms of 
the invariant sections s^^Sy. What then is the set /^ _1 (0) n ^ _1 (0)? 

At a point (xi,X2, • • • ,£12), the condition v = requires that six of the x^s 
are simply negatives of the other six, and then fi = forces the six remaining 
ones, say {x\, X2, x$, X5, xq, x%) , to form a tetrahedron, shown schematically in 
figure 5. Recall that the lengths of the vectors are fixed integers a, b, c, d,e, f . 
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Figure 5: Schematic configuration of vectors 

We have assumed that the numbers a,b, . . . , f satisfy the triangle inequalities 
in triples (otherwise the 6j -symbol is simply zero), so the faces of this triangle 
can exist individually in M 3 . However, it is still quite possible that there is 
no Euclidean tetrahedron r with sides a, b, . . . ,/. The sign of the Cayley 
polynomial V 2 (a? ,b 2 , . . . ,/ 2 ) (whose explicit form is irrelevant here) is the 
remaining piece of information needed to determine whether r is Euclidean, 
flat or Minkowskian. In the last case, we see that /u _1 (0) n ^ _1 (0) = 0, and 
so have proved the second part of the main theorem: that if r is Minkowskian 
then the 6j -symbol is exponentially decaying as k — > oo. 

Suppose on the other hand that V 2 is positive. Then we can find a set of six 
vectors a, b, c, d, e, / in R 3 forming a tetrahedron oriented as shown in figure 6. 
Let r denote both this tetrahedron and the corresponding point 

(a, 6, c, -c, d, e, -e, /, -a, -/, -d, -b) G M 

and t' be its mirror image (negate these 12 vectors). (Of course the whole 
tetrahedron is determined by just three of the vectors, say a, c, e.) It is clear 
that the localisation set ^ _1 (0) n ^ _1 (0) will consist of exactly two G-orbits 
Gt, Gt'. 




Figure 6: Actual configuration of vectors 
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Remark The symbols a, b, c, d, e, f now denote both vectors and their integer 
lengths at the same time! This ought not be too confusing, as it should be clear 
from formulae what each symbol represents. 

Returning to the integral /, since both s k s k are invariant under the diagonal 
action of G, we can apply theorem 5 with respect to the diagonal action of G, 
and obtain an integral over an 18-dimensional manifold Mq: 

1= I (s* i :st)k 12 n = k 9 (t) 3/2 [ (s^sDo-SIg (10) 

where in the right-hand integral, s^Sj, are the descendents of s^,s u , and the 
only thing depending on k is the integrand, which is the kth power of something 
independent of k . The function a is the function on Mq giving the volume of 
the corresponding G -orbit in M , and we view it as part of the measure in the 
integral. 

Let us define ip = log(s fl ,s l/ ) on Mq and tp = log(s^,s u ) on M, so that the 
remaining problem is to compute 

/'= / e k ^o-n G . 

Jm g 

Since the modulus of e k ^ localises on the set /x _1 (0) n ^ _1 (°) = GtUGt' C M, 
the above integrand localises to the two points [r], [t'\. 

5.5 Tangent spaces and stationary phase 

The diagonal G-action does not commute with the other ones, so that Mq 
will not have any kind of induced actions of G 4 or G 6 , but we don't need this 
for the localisation calculation to go through. We always work on the upstairs 
space M not Mq , precisely because the presence of the group actions defining 
the invariant sections being paired is so useful. The tangent space T T M is 
24-dimensional, and contains two 12-dimensional subspaces ker dfi and kerd^ 
which meet in the 3-dimensional space qt. (This degree of transversality can 
be checked explicitly from formulae below, but it should be clear from the fact 
that there are just two isolated critical points in Mq.) Together they span the 
21-dimensional T t Mq, which is orthogonal to iqt. 

Projecting to 18-dimensional T^Mq, we see two 9-dimensional subspaces we 
shall call and W u (the projections of kercfyx and kerdi^) meeting trans- 
versely at the origin. 
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We want to examine the behaviour of tp (its gradient and Hessian) at the point 
[t] G Mq. Let us choose orthonormal bases {Xi, X2, • • • , Xg}, {Yi, Y2, . . . , lg} 
for the transverse 9-dimensional tangent spaces W^^Wy inside T^Mq. Then 
we can need to compute quantities such as Xiip and XjYjif) (also at [r], of 
course!). These can be computed by choosing arbitrary lifts of the vectors to 
T T M and applying them to the G-invariant function ip = \og(s^,s u ) on M. 
This is important, because it is very hard to write down any explicit horizontal 
lifts which would be needed to do computations directly in T^Mq. 

So, to compute something like Xiip one can choose any lift inside ker dfi in 
T T M , and write 

Xiip = Xiip = Xi log(s M , s u ) 

= Sj^dVx.Sn, s u ) + (s M , V x .s v )). 

Now Xi is a generator of the G A action under which s M is invariant, and 
therefore 

which vanishes at r. (Here ri(Xi) really denotes /u(£j) for the Lie algebra 
element £j corresponding to Xi.) For the second term we must first express Xi 
as a linear combination of the Yj and JYj (which span T^Mq), then we can 
lift and use the quantization formula to compute. 

Therefore, introduce the 9x9 matrices Pij and Qij according to 

Xj = y~] PikYk + ^ QikJYk- 

Multiplying by J we get 

j^q = — y~] Qjfc^fe + 5^ PikJYk- 

By applying uJc(Xj, — ) and similar operators to these equations one obtains 

= B G (X t , Yj) = -u G {Xi, Yj). 

These, together with the fact that the bases are orthonormal and span isotropic 
subspaces, determine completely matrices for B and u) on T^(Mc). By a 
similar procedure one can invert the relations: 

Yi = PkiXk - QkiJXk 

jYi = ^2 QkiXk + ^2 Pki J Xk ■ 
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A final point is that since {JQ, JAj} and {Y{, JYi} are both complex-oriented 
orthonormal bases for T^(Mq), the change of basis matrix is special orthogo- 
nal, and hence 

P T P + Q T Q = 1 QP T = PQ T Q T P = P T Q. 

Now we may rewrite the tangent vectors appropriately, lift everything to T T M 
and then apply them to tp via the fundamental formula (recall (— , — ) is con- 
jugate linear in the second factor). For example 

X$ = - 2iti P ikv{Y k ) - 2vr ^ Q*kv(Y k ). 

This right hand side vanishes at f , so indeed Xiip = there. The companion 
formula is 

Yij> = 2mY,P^{X k ) + 2ir^2Q kt Li(X k ) - 2mu{Y). 

Together, these show that is stationary at [r] G Mq, just as in the warm-up 
example. 



5.6 Computation of the Hessian 

Another application of the above formulae, remembering that 

Xfi(Y) = dfi{Y){X) = lo{Y,X) 
will obtain formulae for second derivatives such as 

XjXitp = 2mw G {Xi,Xj) - 2m £ P ik u G (Y k , X,) - 2vr ^ Q ik to G {Y k , Xj) 

where everything in this formula is to be evaluated at r (for example, the first 
term now dies), and we have used the defining identity u>(X,Y) = u> G (X,Y) 
to replace u> by uj g and remove the tildes from the right-hand side. 

We can compute three similar formulae for the second derivatives and form the 
Hessian matrix for ip with respect to the basis {Aj, Yj} of T^Mq : 

(-2m) ( PQT ~ lQQT Q \ 

{ 2m) { Q T P T Q - iQ T Q ) 

We can extract the matrix 

'Q T 
Q 
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form the right, and expand 

, JP-iQ 1 \ 
det [ 1 P T -zQ T ) 

as 

det((P - iQ){P T - iQ T ) - 1) = det(-2QQ T - 2iPQ T ) 

using properties of P and Q discussed earlier. Hence this temporary "unnor- 
malised Hessian" of ip is: 

(-27rz) 18 .(-2i) 9 . det(P - iQ). det(Q) 3 

The reason for separating the last two parts is that det Q is real, whereas P—iQ 
represents the change of basis between {Xi} and {Yj} as bases of T^Mq as 
a 9-dimensional complex vector space (ie, it is the matrix of the hermitian 
form, (P — iQ)ij = HG(Xi,Yj)), so is unitary and contributes just a phase as 
determinant. 

To normalise the Hessian we must compute the volume of the basis {Xi,Yj} 
with respect to the form o~tu G /9\ on T^Mq- Expand, using the shuffle prod- 
uct, the expression (co G /9\)(Xi, X 2 , ■ ■ ■ , Xg, Y±, Y 2) • • ■ ,Yg): since the spaces 
spanned by the Xi and by the Yj are isotropic, the terms appearing are 
simply all possible orderings of all possible products of 9 terms of the form 
uj(Xi,Yj) (the X before the Y). Reordering these cancels the denominator 
9! and we obtain simply det(Q). So the determinant we actually computed 
was (voI(Gt) det(Q)) 2 times what it should have been when computed in a 
unimodular basis. Therefore 

Hess [T] (V0 = (-27ri) 18 .(-2i) 9 .det(P-iQ).det(Q).vol(Gr)- 2 . 

This is the end of the general nonsense. To go any further we have to choose 
explicit bases, although not for T^Mq-, because of the difficulties already men- 
tioned in writing down any vectors there. In the next two sections we will write 
down nice vectors "upstairs" in T T M and show how to lift the computations of 
det(P — iQ),det(Q) into this space. 

5.7 The modulus of the Hessian 

We need to compute det(Q) , where Qij = —L0G(Xi,Yj), and the Xi, Yj are 
orthonormal bases as chosen above. 

Let us start by introducing some useful vectors in T T M , with which to compute 
"upstairs". We make an explicit choice of basis for each of the 12-dimensional 
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spaces ker dfx, ker dv inside T T M . Recalling that they intersect in the space qt, 
we arrange for a suitable basis of this space to be easily obtained from each. 

Let Tl be the infinitesimal rotation about the vector v, acting on the Ith triple 
of vectors from {x\, x 2 , ■ ■ ■ , £12) ■ For example, at any point (x±, X2, ■ ■ ■ , £12) , 
we have 

T] = (bxi 1 ,»x x 2 ,v x x 3 , 0,0, 0,0, 0,0, 0,0,0). 

This vector clearly preserves the condition x± + x 2 + £3 = 0, as well as the other 
three R 3 -coordinate parts of fi. 

Recall that a, c, e are three vectors defining the tetrahedron r . Since a, c, e are 
linearly independent, the vectors T^,T^,T^ span the tangent space to {x\ + 
£2 + £3 = 0} inside the product of the first three spheres of M . Combining four 
such sets of vectors, we see that the 12 vectors T l v , where / = 1, 2, 3, 4 and v is 
one of the three vectors a, c or e, span the space ker d/x at r. For convenience 
these vectors will also be numbered 

rp rp rp rpl rp{ rp{ rpl rpA 

■L 1 ) 2 > • • ■ >^12 J n!-'ci-'e!' l or" ) J e - 

Note that although the formula defines a vector field everywhere on M , we only 
need the tangent vectors at two specific points, namely r and r'. 

We can easily obtain a basis for the infinitesimal diagonal action of G from 
these: 

Ra = T l a +T 2 a +Tl + T* 

is the infinitesimal rotation of all 12 coordinates about a, and similarly we may 
define R c ,R e , each a sum of four 'T's, which together span gr. We will also 
denote these by 

Rl,R 2 , R3 = R a , Rci Re- 
Let u denote an edge of the tetrahedron r, one of the vectors a,b,c,d,e, f . 
Let UJ^j be the infinitesimal rotation about w acting on the pair of spheres 
corresponding to u. For example, if u = a then we have at (xi,x 2 , ■ ■ ■ , x\ 2 ) 

11% = (wx si, 0,0, 0,0, 0,0, 0,io x (-xi), 0,0,0), 

This vector preserves X\ + xg = and hence v, and so do the other U™- We 
want just two vectors wi,w 2 such that U£ 1 ,U® 2 span the tangent space to 
the orbit of G acting on the first and ninth spheres in M at r (compare the 
previous case with the three 'T's.) Projecting into the first and ninth spheres, 
r becomes (a, —a) and U% becomes the tangent vector (w x a,w x (—a)). So 
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all we need to do is pick 101,102 such that a,w\,w 2 are linearly independent. 
In this way we can construct 12 vectors spanning kerdf at r. 

Unfortunately there isn't a totally systematic way of deciding which two values 
of w we should use, given u. We can at least choose them always to be two of 
the three vectors a,c,e, which forces for example the use of U£, U£ among our 
12 vectors (because U£ = 0). 

The twelve explicit choices are as follows: 

i/i, u 2 ,..., u 12 = u a c , u a e , ul u\, u c e , u c a , u d c , ut K, w c , u£, ui 

The same three diagonal generators R a ,R c ,R e can be expressed in terms of 
these vectors by observing that 

R a = K + u b a + u c a + u d a + u e a + ui 

that U% = and that Jj[ (which is the only other not among our chosen 
Ui, Ui, . . . , U12) satisfies u[ = —u£, because the fact that the three sides of 
the tetrahedron r satisfy — e + f — a = implies 

ui + ui = ui - u f f + ui = ui_ f+a = ^ = 0. 

Similarly, one obtains U\ = —U h a and U£ = +Uf}, and hence: 
R a = U 3 + U 6 + U 8 + U 9 -U 12 
Rc = Ui-U z + U 7 + U w + Uu 
R e = U 2 + Ui + U 5 + U 7 + U 12 

In the following calculation, a symbol such as det co ({X i} ; {Yi}) , where {X{\ 
and {Yi} are some sets of vectors, will mean the determinant of the matrix 
whose entries are all evaluations of to on pairs consisting of an element from 
the first set followed by one from the second set (arranged in the obvious way) . 
In the case where the two sets of vectors are both bases of some fixed vector 
space, the symbol det({Xj}/{!^}) will be the determinant of the linear map 
taking Yi 1— > . The grossly- abused subscript i below stands for the complete 
list of such vectors (there are twelve 'T's, three 'i?'s, etc.) We regard all 
vectors as living in T T M , in particular the original orthonormal bases are lifted 
horizontally into it. Let {ei,e 2 ,es} be some orthonormal basis of gr. 

By extending the orthonormal sets of vectors and then changing bases inside 
the spaces ker<i/U and kerdv to bring in the 'T's and 'C/'s , we have 

det(Q) = -det u({Xi};{Yi}) 

= - det uj ({X i, ei,, J ei};{Yi, a, J ei}) 

= - det({T t }/{Xi, ei})- 1 det({U t }/{Yi, ei})- 1 

x detuj({Ti, Jei}; {U, JeJ) 
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The remaining det a; term can be simplified further. Replace three of the T's 
(T 10 , Tn, T 12 ) and three 'f/'s (C/g, C/n, L^) by R\,R 2 ,R 3 . According to the 
earlier expressions for the Ri, each replacement is unimodular, and there is no 
sign picked up in reordering the 'C/'s to put Us,Un,U4 (in that order) last. 
Then we change the 'i?'s back to 'e's and remove them. This gives: 

det co({Ti, Jei}; {Ui, Je^}) 
= detect, ... ,T 9 , Ri, j ei }; {U u U 2 , U 3 , U 5 , U 6 , U 7 , U 9 , U w , U 12 , R t , JeJ) 
= det({^}/{e,}) 2 

x det w({Ti, ... ,T 9 , a, Je^; {U 1 ,U 2 , U 3 , U 5 , U 6 , U 7 , Ug, U 10 , U 12 , e i: Jo}) 
= det({^}/{ ei }) 2 det uj{{T u ... , T 9 }; {U U U 2 , U 3 , U 5 , U 6 , U 7 , U 9 , U w , U 12 }) 

This remaining 9x9 determinant has to be done by explicit calculation. For- 
tunately the good choice of vectors helps enormously. The 'T's have only three 
non-zero coordinates (out of 12), the 'f/'s have only two, and these must overlap 
if there is to be a non-zero matrix entry. So a representative non-zero matrix 
element is something like 

^(Tl,K) = [x.(v x x).(w x x)} = -^[x.v.w] 

where x is whichever of the 'xj's corresponds to the overlap. (It will be plus 
or minus one of a, b, c, d,e, f , depending on whether the overlap of coordinates 
happens in the first or second of the two non-zero slots of the l U '-vector, 
respectively.) 

Writing down the matrix with rows corresponding to T\, T 2 , . . . ,Tq and columns 



1 

47T 



sponding 


to f/ c a , 


jja jjb jjc jjc 

u e ' u a' u ci u a' 


jjd jje 

u c i u a ' 


u e c ,ui 


gives 








( ° 





o 


cae] 

















\ 





[ace] 


[bca] 




















[aec] 





[6ea] 





[cea] 

























[cae 





[dac 





[eac] 

























[eca] 






















— [cea 


] [dec] 


































— [eac] 


[fae 







— [ace 














— [eca] 





[fee 




\— [aec] 


























/ 



where [ace] is just a shorthand for the vector triple product [a.c.e]. Substituting 
the relations b = —a — c,d = c — e, / = a + e and extracting the factor of [ace] 
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gives 



[ace] 
Air 











-1 



















1 























-1 





-1 





1 























1 





-1 





1 























-1 




















-1 



































-1 








-1 














1 





1 


Vi 























0/ 



This matrix has determinant [ace] 9 / (An) 9 . 

The various change-of-basis determinants may be evaluated in terms of orbit 
volumes. If we denote by sgn({Aj}/{Yj}) the sign of the determinant of the 
appropriate transformation then 

det({T t }/{Xi, a}) = S ga{{Ti}/{Xi, e t }) vol^Tj 

because {Xi,e{} is £? r -orthonormal. The volume is given by 

vol(G 4 r) r „ 



vo\ BT {Ti} 



voL(G) 4 



ace 



by lemma 6 and the fact that the 12 T's separate into four orthogonal triplets 
coming from the Lie algebra elements a, c, e G M 3 . Similarly we have for the 
'{/' case: 

det({U t }/{Yi,e t }) = sgadUij/iY^e^volB^Ui} 

The volume term may be expressed via lemma 7. This will involve a product 
of six terms of the form vol p {£j}: for each choice of u, we have to calculate 
the area spanned by the vectors w\ and wi once projected into the orthogonal 
complement of u. This is just | [w\.W2-u]/u\ . Substituting the explicit choices 
we made gives 

vol(G 6 r) : 



vo\ Br {Ui} 



ace 



vol p (G/T)6 Y[a 

where the product on the right hand is simply abcdef . Yet another application 
of lemma 6 yields: 



det({^}/{ej) 5 



vol(Gr) 
vol p (G) 



ace 



The sign terms here depend on the original choice of orthonormal bases, which 
we did not specify. So we do not yet know the actual sign of det(Q). Similar 
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terms will appear in the computation of det(P — iQ) , however, so that the 
product of the two terms does not depend on the original choice. So far we 
have 



det(Q) = -sgn({T i }/{X i ,e i }) S gn({U i }/{Y i ,e i }) 

/ vol(G 4 r) 14VV vol(g 6 r) [ace^ y 1 / vol(Gr)[ace] \ 2 
X \^(W [aCel ) Ul p (G/T)6 n« / V vol p (G) J 

(11) 



5.8 Phase of the Hessian 



We work out det(P — iQ) (defined by ((P — iQ)ij = Hc(Xi, Yj)) using similar 
techniques. We choose slightly different bases in T T M this time. 

For each face of the tetrahedron r, numbered by I as earlier, choose three 
infinitesimal rotation vectors T l v by letting v be an exterior unit normal vector 
vi to the face or one of the two edges xsi-2,%3i-i of that face. These occur in 
clockwise order, so that 



with A\ twice the area of the Ith face. (Of course at the point r , we know that 
each Xi is just one of the vectors a, b, c, d, e, f or their negatives. However, it is 
easier to calculate without substituting these yet.) Pick a set of 12 vectors 
rather as before, except that given an edge u, we allow w to be the exterior 
unit normal to either of the two faces incident at u . Order the two choices so 
that the first cross the second points along u (in fact this corresponds to V{ 
coming before Vj iff i < j). Figure 7 shows these where these vectors are in 
R 3 , given the tetrahedron r. We will refer to the chosen vectors as U[, . . . , U[ 2 



det(P - iQ) = det HodYi}; {X,}) 

= det H({Yi, a}; {X u a}) 

= det({^}/{^, e,})- 1 det({Ui}/{Y h e,})" 1 det H{{U'^ {T>}). 



Since we are know the determinant is actually just a phase, we can throw away 
any positive real factors appearing during the computation. For example, the 
correcting determinants above may immediately be replace by correcting signs, 




XZI-2 x Xzi-l = Am 
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because the difference (a volume) is positive. This principle also facilitates the 
direct computation of det i?({T^}; {U[}) too. 

Let us compute sample non-zero elements (once again, most of the matrix 
elements will be zero): 

H{Ul,T l v )=u{TlJUZ + iUZ) 

1 i 

[x.(v x x).(x x (w x x))] + - — k[x.(v x x).(w x x)\ 



4irx 3 4nx 2 
= —(x(v.w) + i[xvw]) 

using the earlier notation for the triple product, and with x being whichever 
of the 'xj's corresponds to the overlap of the non-zero coordinates of £7^,7^. 
There has been some simplification because w.x = 0. 

Let us immediately forget about the 4-7T factors. If the ith and jth faces meet 
in an edge u , oriented along the direction of Vi x Vj , then the exterior dihedral 
angle, written 6 U or 6ij , is defined (in the range (0,7r)) by 

usin(6ij) = [u.Vi.Vj] 
cos(9ij) = Vi.Vj. 

In performing the computation we run across three kinds of non-zero matrix 
elements. If v = v\ is normal to the lib. face, and w = then we obtain 
= |x|e* efci . (If k = I it is simply |x|.) If instead v is a vector lying in 
the Ith. face, and w = v^, then we distinguish according to whether k = I or 
not. In the case of equality, we get case we get 0, iAi or — iA\ according to 
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whether v is x, its successor, or predecessor in the anticlockwise cyclic ordering 
around the face. If k 7^ I then we get 0, iAia^i or iAia^i, according to the 
same conditions. 



We can throw out the area factors and the eight powers of i coming from these 
second and third cases, and end up with a matrix: 



/ a 


aai3 


b 


ban 


c 


CQ12 

















\ 








-1 


— ai4 


1 


Q12 




















1 


ai3 








-1 


— ai2 
































CQ12 


c 


d 


da24 


e 


e«23 


























-1 


— OI24 


1 


Q23 




















a\2 


1 








-1 


— Q23 








actri 


a 




















e«23 


e 


/ 


f a 34 


Ql3 


1 


























-1 


—034 


-a 13 


-1 




















e*23 


1 














ban 


b 








da-2.4 


d 








/Q34 


/ 








e*i4 


1 








— «24 


-1 














V 





— ai4 


-1 




















«34 


1 J 



Easy row operations and then permutation of the rows reduces the determinant 
to minus that of the direct sum of the six 2 x 2 blocks 



Discarding the sines, which are positive and so only affect the modulus of the 
determinant, we find that its phase is 



Let us collect up det(P — iQ) det(Q) finally. The annoying sign terms may be 
combined into 



The first term is a product of four signs arising from orientations of a three- 
dimensional vector space, and the second a product of six arising from two- 
dimensional spaces. All these signs are positive (an easy check). 




The determinant of such a block is 



-a.ij.1i sin(0jj). 



sgn({T,}/{7;'})sgn({^}/{^})- 
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Thus we have for the Hessian: 

9 



Hess [r] (V) = -i(27r) 18 2 9 e^ ea 



[ace\ 
4tt 



vol(G 4 r) r l4 ^ 1 ( vol(G 6 r) \acef\ 1 { \acJ ^ 2 



vol p (G) 4L V Vvol p (G/T)6 [la 7 Vvolp(G) 



(12) 



Looking at the argument again, it is easy to see that the Hessian at [r'] is the 
complex conjugate of this one. 



5.9 The overall phase of the integrand 

We need to account for the 0-order contributions Y>([t']), ^([ r D of the integrand 
at the two critical points. Since the two sections being paired were fixed to 
have norm 1 along their critical regions, these 0-order contributions also have 
modulus 1. We start by calculating the phase difference. 

The chosen lifts r, r' of these points lie on the slice \x = inside M . Since this 
slice is just a product of four "spaces of triangles", each of which is a principal 
G-space, (G = 50(3)) there is a unique element of G A which translates r to 
t' . In fact it is easy to describe such an element g = (g±, 172, #3, 54) explicitly. 
The element g\ must rotate the triangle (a, b, c) (the projection of r into the 
first three sphere factors of M) to its negative (— a, — b, — c). Therefore it is 
the rotation of tt about the normal to the triangle's plane. The other gi are 
similarly half-turns normal to their respective triangular faces. 

We must compare the values of the pointwise pairing L, s u at r, r' . Define for 
each face a lift gi of gi into SU (2) by lifting the path of anticlockwise rotations 
from to tt. Together these form g 6 SU{2) A . Since s M is S*C/(2) 4 -invariant: 

s^(t') = s^gr) = gg^is^gr)) = 5((5^ 1 s M )(r)) = <?(s M (r)) 

By contrast, s u is not SU(2) 4 -invariant, though it is SU(2) 6 -invariant. We 
can write an equation like the above but we need to know what (g _1 Sj,) is to 
perform the last step. Now s u is a sextuple tensor product, and we can study 
the action of g^ 1 on it by looking at the action on the six factors individually: 

F% = (~gi 1 ,~93 1 )s aa ® (gi\gl l )s bb ® ■ ■ ■ ® (fifc^rV 
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Using the diagonal invariance of each section of the form s aa , we have identities 
like 

Now <7i,<?3 are lifts of rotations through tt about directions normal to the two 
faces of the tetrahedron meeting at side a, so the composite g^gi is a lift of 
an anticlockwise rotation through an angle equal to twice the exterior dihedral 
angle 9 a about the vector a, but it is slightly tricky to decide which lift. Let 
us denote by r a the lift of the path of anticlockwise rotations from to 29 a , 
and write g% l gi = S a f a , for some S a = ±1. Then 

((sr\23~V a )(a, "«) = (lAraX^KW^Xa, -a))) 
= (l,<5 a f a )( S a >, -a)). 

The action of r a on the fibre of the bundle C® a —> P 1 at —a is multiplication 
by e~ m9a (remember that it acts as e l9a on the fibre of the tangent bundle C® 2 
at a). The sign 5 a acts as its ath power 5°. So, using the invariance of the 
hermitian form, 

{s»{t%~s v {t')) = (~9(~s,(r)),~9((U S a) e^^s v (r))) 

= {U 6 a) e^^(^(r), S ~,(r)). 

In fact this identity is independent of the choice of lifts cji . For example, chang- 
ing the lift of g± negates 5 a ,Sb,5 c , changing the right-hand side by (— l) a+6+c , 
which is +1 because of the parity condition on a,b,c. Further, if we imag- 
ine varying the dihedral angles of the tetrahedron in the range (0, 7r), all our 
chosen lifts are continuous and so we may evaluate the sign (FJ <5") by defor- 
mation to a flat one, for example where 9 a ,9b,9 c are ir and the others 0. In 
this case, 52 = 53 = 54 = ~9\ , and so all the S a turn out positive. Hence 
VKb"']) = e lk ^ a9a i{j([T]). Because the 6j-symbol itself is real, the two values 
of ip must be conjug ate. Therefore = ±ei ik ^ a0a , but we still have an 

annoying sign ambiguity. 

There are really three separate sign problems: how the sign depends on k (for 
fixed a,b, . . . , / ); how it varies as we alter a,b, . . . , /; and one overall choice 
of sign. One step is easy: the phase conventions on the sections used in the 
pairing implied for example that s ka ' ka = ( s aa )® fc and that the integrand was a 
kth power of another function, the sign above must actually be (— l) k or or 1. 

To do better requires a frustrating amount of work, which will only be sketched 
here. Recall that the signs in definition 3 were fixed by writing down explicit 
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polynomial representatives of the trilinear and bilinear invariants. If we remove 
a suitable branch cut from each copy of the sphere in M , leaving a contractible 
manifold, we can extend this definition to allow real values of the variables 
a,b, . . . , /, and extend VKM) to a real-analytic function of these variables (at 
least locally). It is obtained by pairing two holomorphic sections of a trivial 
bundle with a hermitian structure which still satisfies the quantization formula 
(4). There is another way of computing the phase difference above, based on 
choosing two paths r — > r' in M, one on which fi = and one on which 
v = 0, and computing the holonomy around the resulting loop. This can be 
done by computing the symplectic area of a bounded disc. To carry this out 
appropriately one must be very careful with which disc: once the symplectic 
form no longer has integral periods on H2(M) , this matters. Further, the most 
obvious paths and disc intersect the branch cuts, so one must account for this 
too. Ultimately one obtains an analytic expression 

Y>([t']) = ±e^^* + '^ OT 

where the analyticity restricts this sign to a single overall ambiguity. (Note that 
at integral values of the lengths, we can see the sign (— l)^ a appearing.) One 
could compute this sign using a single example, but as the reader will judge 
from this terse paragraph, the author is so bored with fixing signs that he no 
longer cares to! The experimental evidence in [14] confirms that the sign is 
positive. 

Therefore 

y,([r'] = (_i)2> e 3<*2:«0« and V(M) = (-1) E a e~i ik ^ ada . (13) 



5.10 Putting it all together 

We combine the original integral definition (7) with the asymptotic normali- 
sation factors (8), (9), the reduction (10), the stationary phase evaluation (6) 
incorporating the Hessian (12) and 0-order terms (13). 



ka kb kc 
kal ke kf 



(j[Vka~+tj 



2tt 
~k 



k 
2 



(vol(^ 1 (0))vol(i/- 1 (0)))-U 9 



k\ 2 



|lfcyj 0,0a 



+ 



e2 



Hessr Tl (^) 



Hessr T /i (ip) 
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The terms from the Hessian and the normalisation involving 

vol(/i -1 (0)) = vol(G 4 r) = vol(GV) and vol(i/ -1 (0)) = vol(G 6 r) = vol(GV) 

cancel. The normalisation factor (n(^ a + I)) 5 cancels with the term in the 
Hessian involving Y\ a > contributing asymptotically simply A; 3 . What remains 
is 

(2^)t2^fc~f [ace]~^ vol p (G)" 1 vol p (G/T)- 3 cos|^(A;a + l)y + ^ j. 
Substituting in the volumes 8ir 2 and Air of G and G/T gives 



ka kb kc 
kd ke kf 



v3? cos {5>+i>t+j} 



where V = |[ace] is the (scaling-independent) volume of r. This completes 
the proof of the theorem. 



6 Further geometrical remarks 

6.1 Comparison with the Ponzano— Regge formula 

It is important to note that the formula (2) is not the same as the original 
Ponzano-Regge formula. There are two main differences, apart from the trivial 
fact that they label their representations by half-integers instead of integers. 

Their claim, in our integer-labelling notation, is that for large a, b, c, d,e, f: 



a b c 
d e f 



t^cos jj^(a + l)^ + ^| if t' is Euclidean, 
exponentially decaying if t' is Minkowskian, 

(14) 



where r' is a tetrahedron whose edges are a + 1, 6+1, . . . , / + 1 and whose dihe- 
dral angles 0' a and volume V are therefore slightly different from those of our 
r. This difference is worrying, as it is quite possible to find sextuples of inte- 
gers such that t' is Euclidean yet r is Minkowskian, in which case the formulae 
seem to conflict: is the 6j'-symbol exponentially or polynomially decaying in 
this case? 



Geometry & Topology, Volume 3 (1999) 



62 



Justin Roberts 



The second difference explains this. The Ponzano-Regge formula (14) is only 
claimed as an approximation for large irreducibles, rather than an asymptotic 
expansion in a strict sense as in theorem 1. Therefore the only meaningful 
comparison between the two formulae is to examine how their function behaves 
as we rescale a, b, c,d,e,f by k — > oo in the precise sense of our theorem. 

Although for small k it is possible that r' might be Euclidean when r is not, 
eventually the shift in edge-lengths becomes insignificant and either both are 
Euclidean or neither is. Therefore there is no inconsistency between cases in 
the two formulations. 

As for comparing the actual formulae, the asymptotic behaviour of the Ponzano- 
Regge function is 

because V and V agree to leading order in k . The only problem is the dihedral 
angles relating to slightly different tetrahedra. Fortunately we may easily show 
that 

e ikZ(a+i)0' a ^ e ikY.{a+i)e a 

by applying the Schlafli identity (see Milnor [11]), which says that the differ- 
ential form adOa vanishes identically on the space of Euclidean tetrahedra. 
Therefore there is no inconsistency. 

Remark The case of a flat tetrahedron is not covered by either formula. 



6.2 Regge symmetry and scissors congruence 

Suppose one picks out a pair of opposite sides of the tetrahedron denoting the 
6j -symbol (as in figure 1), say a,d. Let s be half the sum of the other four 
labels (twice their average). Define: 



a' = a 



b' = s-b 
d' = d d = s — c 



(16) 

e = s — e 



f' = s-f 
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Regge discovered that the 6j -symbols are invariant under this algebraic op- 
eration (the easiest way to see this is to look at the generating function for 
6j -symbols, [15]): 



We can also consider this as a geometric operation on a tetrahedron, altering its 
side lengths according to the above scheme. It is not meant to be obvious that 
the result of applying this to a Euclidean tetrahedron will return a Euclidean 
one! 

Regge and Ponzano considered the effect of this symmetry on the geometrical 
quantities occurring in their asymptotic formula, mainly as another check on 
its plausibility. They discovered that the volume and phase term associated to 
a Euclidean tetrahedron are indeed exactly invariant. This is amazing, given 
that it would be consistent with their appearance in an asymptotic expansion 
for them to change, but by lower-order contributions. 

Let us reconsider this surprising geometric symmetry. First note that the sym- 
metry is an involution: if one thinks geometrically, it corresponds to reflecting 
the lengths of the four chosen sides about their common average. These involu- 
tions, together with the tetrahedral symmetries, form a group of 144 symmetries 
of the 6j-symbol, isomorphic to S4 x S3 (see [15]). 

V. G. Turaev pointed out to me that the term ^ liOi in the phase part of 
the formula (2) for a tetrahedron r is reminiscent of the Dehn invariant <5(r). 
Actually it would be fairer to say that it is the "Hadwiger measure" (or "Steiner 
measure" ) fi\ (r) . Both invariants are connected with problems of equidissection 
of three-dimensional polyhedra. 

Two polyhedra are scissors congruent if one may be dissected into finitely-many 
subpolyhedra which may be reassembled to form the other. (Hilbert's third 
problem was to determine whether three-dimensional polyhedra with equal vol- 
umes were, as is the case in two dimensions, scissors-congruent. Dehn used his 
invariant to solve this problem in the negative.) 

The modern way of looking at the problem is to define a Grothendieck group 
of polyhedra V . We take Z-linear combinations of polyhedra in R 3 with the 
relations: 




PuQ = P + Q-Pf]Q 

P = if P is degenerate 

P = Q if P, Q are congruent 



(17) 
(18) 
(19) 
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Volume is an obvious homomorphism V — ► M. The Dehn invariant is a less 
obvious one V — > R <8>z (M/7rZ) , denned for a polyhedron by summing, over its 
edges, their lengths tensor dihedral angles: 



Sydler proved that these two invariants suffice to classify polyhedra up to 
scissors-congruence: two such are scissors-congruent if and only if they have 
the same volume and Dehn invariant. See Cartier [2] for more details. 

If we look for homomorphisms V — ► R which are continuous under small per- 
turbations of vertices of a polyhedron, then volume is the only one (up to 
scaling). However, if we remove the condition (18) on degenerate polyhedra 
from the axioms defining V, there is a four-dimensional vector space of con- 
tinuous homomorphisms, spanned by the following Hadwiger measures (picked 
out as eigenvectors under dilation): 



See Milnor [10] or Klain and Rota [8] for more on these beautiful functions. 
The relationship with the Regge symmetry of tetrahedra is as follows: 

Theorem 9 The Regge symmetry (16) takes Euclidean tetrahedra to Eu- 
clidean tetrahedra, preserving volume, Dehn invariant and Hadwiger measure 
ll\ . (Remark: simple examples show that Regge symmetry does not preserve 
the surface area measure ji2-) 

Proof The tour-de-brute-force of trigonometry in appendices B and D of [14] 
contains all the calculations necessary to prove this. They demonstrate that 
under the Regge symmetry, which is a rational linear transformation A of 
the six edge lengths, the dihedral angles also transform according to A. The 
orthogonality of this matrix and the fact that we may view the Dehn invariant 
as being in R (g)Q (R/ttZ) = M ®z (R/ttZ) demonstrate its invariance, as well 
as that of \i\ . The volume is checked by straightforward calculation using the 
Cay ley determinant. □ 

Corollary 10 The orbit under the group of 144 symmetries of a generic tetra- 
hedron consists of twelve distinct congruence classes of tetrahedra, all of which 
are scissors-congruent to one another. 



M3(P) 

MP) 
MP) 

MP) 



vol(P) 
lvol(dP) 

x{P) (the Euler characteristic) 



(20) 
(21) 
(22) 
(23) 
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Remark The fact that a tetrahedron is scissors-congruent to its mirror-image 
was proved by Gerling in 1844 (see Neumann [13]), using perpendicular barycen- 
tric subdivision about the circumcentre. One would expect that for the Regge 
symmetry, which is also a "generic" scissors congruence (as opposed to a "ran- 
dom" coincidence of volume and Dehn invariant for two specific tetrahedra), a 
similar general construction might be given. What is it? 

6.3 Further questions 

1. Ponzano and Regge give an explicit formula for the exponential decay of the 
6j -symbol, in the case when no Euclidean tetrahedron exists. Amazingly, it is 
an analytic continuation of the main formula, incorporating the volume of the 
Minkowskian tetrahedron which exists instead, and with the oscillatory phase 
term converted into a decaying hyperbolic function. Can this be extracted from 
a similar procedure? 

2. Can similar geometrically- meaningful formulae be obtained for general spin 
networks, the so-called 3n//-symbols? 

3. The calculation in this paper is comparatively crude, since it computes a 
pairing of 12-linear invariants when one could really do with a pairing of 4-linear 
invariants (see the remark of section 2.2). The space whose quantization gives 
quadrilinear invariants is 2-dimensional, in fact a sphere with a non-standard 
Kahler structure. Can one work directly on this space instead? (Possibly 
any advantage in dimensional reduction is lost when one needs to do explicit 
calculations, which end up like the ones here). 

4. Can similar formulae be obtained for other groups, and does their associ- 
ated geometry have any physical meaning? The 6j-symbols are scalars only 
for multiplicity-free groups such as SU(2). In general they live in the tensor 
product of four trilinear invariant spaces, in which one would need preferred 
bases. 

5. Can one obtain similar formulae for quantum 6j-symbols, which arise as 
pairings in the quantization of moduli spaces of flat connections on the 4- 
punctured sphere? 



Geometry & Topology, Volume 3 (1999) 



66 



Justin Roberts 



References 



R Carter, G Segal, I MacDonald, Lectures on Lie groups and Lie algebras, 
LMS Student Texts 32, Cambridge University Press (1995) 

P Cartier, Decomposition des polyedres: le point sure le troisieme probleme de 
Hilbert, Asterisque 133-134 (1986) 261-288 

W Fulton, J Harris, Representation theory, Springer GTM 129 (1991) 

V Guillemin, S Sternberg, Geometric quantization and multiplicities of group 
representations, Inventiones Math. 67 (1982) 515-538 

N Hitchin, Metrics on moduli spaces, Contcmp. Math. 58 part 1 (1986) 157- 
178 

L Kauffman, S Lins, Temperley-Lieb recoupling theory and invariants of 3- 
manifolds, Annals of Maths Studies 134, Princeton University Press (1994) 

A A Kirillov, Geometric quantization, from: "Dynamical systems IV" , Ency- 
clopaedia of Mathematical Sciences 4 (V I Arnol'd and S P Novikov, editors) 
Springer (1990) 

D A Klain, G-C Rota, Introduction to geometric probability, Cambridge Uni- 
versity Press (1997) 

D McDuff, D Salamon, Introduction to symplectic topology, Oxford Univer- 
sity Press (1995) 

J Milnor, Euler characteristic and finitely additive Steiner measures, from: 
"Collected papers vol. 1", Publish or Perish (1993) 

J Milnor, The Schldfli differential equality, from: "Collected papers vol. 1", 
Publish or Perish (1993) 

D Mumford, J Fogerty, F C Kirwan, Geometric invariant theory, Ergeb- 
nisse der Mathematik 34, Springer (1994) 

WD Neumann, Hilbert's 3rd problem and invariants of 3-manifolds, from: 
"The Epstein birthday schrift" , (Igor Rivin, Colin Rourke, Caroline Series, ed- 
itors) Geometry and Topology Mongraphs, 1 (1998) 383-411 

G Ponzano and T Regge, Semi-classical limit of Racah coefficients, from: 
"Spectroscopic and group theoretical methods in physics" (F Bloch, editor) 
North-Holland (1968) 

D A Varshalovich, A N Moskalev, V K Khersonskii, Quantum theory of 
angular momentum: irreducible tensors, spherical harmonics, vector coupling 
coefficients, 3nj -symbols, World Scientific (1988) 

N Y Vilenkin, A U Klimyk, Representation of Lie groups, and special func- 
tions, vol. 1, Kluwcr (1991) 

E Wigner, Group theory, Academic Press (1959) 



Geometry & Topology, Volume 3 (1999) 



